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Abstract 

The presynaptic free Ca^"*" dynamics of neurons due to various stimulation protocols 
is investigated in a mathematical model. The model includes Ca^"*" influx through 
volt age- dependent Ca^"*" channels, Ca^"*" buffering by endogenous and exogenous 
buffers as well as Ca^"*" efffux through ATP-driven plasma membrane Ca^+ pumps 
and Na"'"/Ca^"'" exchangers. We want to support a specific way of modeling which 
starts on the level of single proteins. Each protein is represented by characteris- 
tics that are determined by single protein experiments and that are considered to be 
widely independent of neuron types. This allows the applications of the model to dif- 
ferent classes of neurons and experiments. The procedure is demonstrated for single 
boutons of pyramidal neurons of the rat neocortex. The corresponding fluorescence 
measurements of Koester & Sakmann (2000, J. Physiol., 529, 625) are quantitatively 
reproduced. The model enables us to reconstruct the free Ca^^ dynamics in neurons 
as it would have been without fluorescence indicators starting from the fluorescence 
data. We discuss the different Ca^+ responses and find that during high-frequency 
stimulation an accumulation of free Ca^+ occurs above some threshold stimulation 
frequency. The threshold frequency depends on the amount of fluorescence indicator 
used in the experiment. 
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1 Introduction 



Ca^+ is an important cellular messenger. In particular, Ca^+ plays a key role 
in synaptic transmission. For example, it triggers the secretion of neurotrans- 
mitter [1,2] and is involved in synaptic plasticity [3] being most likely related 
to the neural bases of memory [4]. Important phenomena in this context are 
presynaptic facilitation, short-term potentiation (STP), long-term potentia- 
tion (LTP) and long-term depression (LTD). It has been shown that the in- 
duction of long-term effects (LTE) are related to a prior Ca^+ influx into the 
neuron [5]. 

There exist several theoretical and numerical approaches to describe intracel- 
lular calcium dynamics. Some models investigate presynaptic Ca^"*" dynam- 
ics in response to single action potentials [6,7,8]. The role of free calcium in 
transmitter release and in synaptic facilitation has been discussed in detail 
[6,9,10,11]. These models focus on the geometry of the presynaptic compart- 
ment and the distribution of calcium transporting proteins in the membrane. 
The role of intracellular organelles with the capability to store and release free 
calcium (e.g. endoplasmatic and sarcoplasmatic reticula) has been analyzed 
[12,13,14]. Also the role of slow and fast intracellular buffers has been inves- 
tigated [15,16]. Indeed, most of the incoming calcium ions are bound to some 
intracellular buffers [17], such that this issue is very important. 

The philosophy of the above models is to focus on the spatial distribution of 
calcium ions in different compartments of the cell. The connection of calcium 
and transmitter release has been investigated mostly in two giant synapses: 
the calyx of Held [18] and the squid giant synapse [19]. For example it has 
been discussed for the squid giant synapse that the calcium enhancement in 
the direct neighborhood of the channels is related to transmitter release and 
facilitation [20]. To this end the channels enter this model as discrete entities. 
We would like to reconsider this approach on the basis of new data that are 
available today. 

Single channel properties have become more and more available [21]. In the 
time of the work of Fogelson and Zucker the exact single protein currents 
and open dynamics were not available. Therefore, these data were modeled 
with some free parameters. The question we would like to address is whether 
we have sufficient information on single protein properties at hand in order 
to explicitly and quantitatively describe the activation and activity of each 
protein involved in the transmembrane calcium flux. 

To this end we will proceed as follows: 

• Instead of describing large squid giant synapse we restrict ourselves in a first 
step to rather small synapses. We focus on the average calcium transient 
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in the synapse and neglect the effect of calcium diffusion. This prevents us 
from considering LTE or synaptic facilitation in detail due to the lack of 
detailed spatial information on inhomogeneous calcium distributions. The 
average calcium transient is nevertheless an important indicator since the 
exact amplitude and time course of calcium, i.e. the exact shape of the 
calcium transient, is believed to be important for transmitter release [22]. 
The average calcium transient is also suitable for the analysis of the quality 
of a model based on single protein properties. 

• The calcium dynamics is described by well known ordinary differential equa- 
tions. Within these equations the single protein properties are isolated and 
treated separately. 

• The model parameters concerning the behaviour of single proteins are de- 
termined on the basis of independent single protein experiments. This set 
of parameters will be called universal because it will be assumed that the 
properties of the proteins do not primarily depend on the neuron under 
consideration. This is surely an approximation which has to be further val- 
idated. 

• The model parameters concerning the specific setup of an experiment, i.e. 
extension of the synaptic compartment, calcium concentrations, stimula- 
tion protocol, use of fluorescence indicator, etc. are adjusted to a specific 
experiment measuring calcium transients in response to various stimulation 
protocols in single boutons of pyramidal neurons of the rat neocortex [23]. 
These specific parameters have to be adjusted for each experiment sepa- 
rately. 

• The remaining parameters (in our case this will be the protein densities in 
the membrane of the neuron) are used in order to fit the data. It is important 
that the number of fit parameters is minimized in order to get a suitable 
test of the model. 

Note, that the separation of universal and specific parameters in the model 
opens the possibility to use the universal parameters, which have been deter- 
mined once in independent experiments, for different experimental setups. On 
one hand, this considerably reduces the number of fit parameters and strongly 
restricts possible outcomes of the model. On the other hand, this allows to 
compare the results of different experiments using different protocols or cal- 
cium indicators. Indeed, it has been mentioned before, that the amount and 
type of calcium indicator changes the calcium transients in a non-negligible 
way [15,16]. Therefore, we consider this model approach as a first step towards 
a link between different experiments opening the possibility of a comparative 
approach to calcium data. 

Our main interest is a quantitative description of presynaptic free Ca^+ dy- 
namics in response to single action potentials and to different stimulation 
protocols. We apply the model to single boutons of the pyramidal neurons in 
the rat neocortex and compare the results to the corresponding experiment 
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[23]. This experiment provides most of the necessary informations in order to 
determine the specific parameters of the model. We will look for characteris- 
tics in Ca^"*" time courses which are specific for the used stimulation protocols. 
Furthermore, we will analyze the influence of the calcium indicator on the re- 
sults of the experiment and discuss the possibility to reconstruct the calcium 
transient as it would have been without indicator. 



2 Methods 



The intracellular free Ca^"*" dynamics due to a predefined, time dependent 
membrane potential of a presynaptic neural domain is described by using a 
set of ordinary differential equations. This neural domain is represented by a 
cell membrane which separates the extracellular space from the intracellular 
space. Specific membrane proteins are embedded into this membrane. Some of 
these membrane proteins and in addition some intracellular proteins determine 
the dynamics of the free Ca^"*" ions concentration [24,25,26]. 

The model geometry consists of a homogenous reaction volume (one-compart- 
ment) which is situated in a homogeneous extracellular space. We are using a 
spherical symmetry of the reaction volume. However, the results remain valid 
for other shapes provided that the surface to volume ratio is conserved. In 
general presynaptic boutons are relatively small, i.e. having a diameter of the 
order of microns [23]. In this study we assume that the process of diffusion of 
Ca^"*" ions is fast enough to reach an equilibrium distribution on a time scale 
of one ms. This assumption has to be reviewed considering larger reaction 
volumes as found in postsynapses. Consequently, also the surface densities of 
membrane proteins have to be interpreted as average values over the whole 
reaction surface. 

Starting from a cell in equilibrium the application of a depolarizing membrane 
potential V induces a free Ca^+ ion current between the reaction volume and 
the extracellular space through diverse membrane proteins. This leads to a 
change of the intracellular free Ca^"^ concentration c = Csh"^ which is rep- 
resented by an ordinary differential equation 

= ^ { Ji (l^ (t) , V (C)) - Je (C) +L} + 5en (c) + ^ex (c) (1) 



in analogy to other studies (see for example [27] eq.(l) or [13] eq.(9)). Every 
source and every sink of the free Ca^"^ ions is determined by its own term. 
Ji is the influx current density per membrane surface unit through voltage- 
dependent Ca^+ channels (VDCC) [28,29,30]. It directly depends on the ap- 
plied membrane potential V (the time dependence in Eq. (1) is noted explicitly 
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to emphasize that V (t) is treated as input signal), and the effective reversal 
potential V and is described in detail in section 2.1. Jc is the efflux current 
density per membrane surface unit caused by ATP-driven plasma membrane 
Ca2+ pumps (PMCA) [31,32,33] and Na+/Ca2+ exchangers (NCX) [34,35,36]. 
It depends on the free intracellular Ca^"^ concentration c and is described in 
detail in section 2.3. The leakage surface current density L is determined by 
the steady state conditions and represents all currents which are not described 
in our model in an explicit manner and that determine the equilibrium state 
(section 2.6). The terms -Ben and -Bex stand for the action of endogenous and 
the exogenous buffers [37,38,39], respectively, which bind to and disengage 
from the free Ca^"*" ions (see section 2.2). The geometry factor G = 0/V is 
the surface to volume ratio. It translates the current densities (Ji, Jc, L ) 
into changes of free intracellular Ca^"'' concentration, z and F are the valence 
charge of Ca^"*" ions and the Faraday constant, respectively. 

Note, that the separation of universal and specific parameters is not reflected 
in Eq. (1). However, in the following sections the explicit form of the currents 
entering the equation will be developed on the level of single membrane pro- 
teins. It is at this point that the parameter classification will be mirrored in 
the mathematical prescription. 



2.1 Ca?+ influx through VDCC 



The VDCC is characterized by a current-voltage relation based on both a 
voltage dependent current through the open pore jsmgic channel ^ 
voltage dependent opening probability [28]. We use the following description 
of the Ca^"*" infiux current density J through VDCC [40] 

Ji (V (t) , V (c)) = g, {V (t)) g, (v (c) - V (t)) . 

^ ^ ^ . ^ (2) 

rsingle channel 
open 



Ji is weighted with the surface density of the channels py being the most 
important cell-type specific parameter. Practically all other parameters defin- 
ing the single channel properties belong to the class of universal parameters. 
The Ca^^ current through the open pore jsj^^gie channel jg driven by the elec- 
trochemical gradient over the membrane [41,42]. A current caused by changes 
of the membrane potential is followed by changes of the intracellular Ca^"*" 
concentration and these retroact on the current. The Ca^+ current due to 
the membrane potential gradient is approximated by a linear voltage-current 
relation with open pore conductivity y^. This is justified for physiologically 
relevant membrane potentials for which the VDCC current-voltage relation 
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has been found to be indeed nearly linear [43,44,45]. The important potential 
difference entering Eq. (2) is the one relative to the Ca^+ reversal potential 
V which dynamically incorporates the Ca^+ concentration gradient over the 
membrane into the model using Nernst equation [46] 

V{c)^^ln{^)-AV.,. (3) 



R is the molar gas constant and T is the absolute temperature. Cext stands for 
the external Ca^+ concentration. AVes corrects the exact reversal potential 
for the linear approximation. As it has been found that efflux Ca^+ currents 
through VDCC at voltages above the reversal potential are negligible [47,48] 
we set /qp^JJ*^ = for V > V. Independently from the definition of the 
reversal potential by the Nernst equation it should be mentioned that its Ca^+ 
dependency has only minor effects on the Ca^"*" dynamics. Note that the exact 
value of AVes may depend on the ion concentrations in each cell. 

The time dependence of the single channel open probability gy is modeled by 
a single exponential approximation 

= (^v (V) - g.) \ (4) 

which may be related to experiment using data for the average number of 
open channels in an ensemble of channels [43,45,49]. The open probability ap- 
proaches its asymptotic value ^v(^) with a time constant r. g^iy) is different 
for each membrane potential and is described by a sigmoidal function 

exp ((14 - ^) y + 1 



being the best approximation to experimental data [43,44,45]. Here, Vh is the 
half activation voltage and k describes the steepness of the asymptotic opening 
probability g^. 



2.2 Ca^"*" buffering via endogenous and exogenous buffers 



The main part of the intracellular Ca^+ is bound to an endogenous buffer 
[17,50,51]. In addition, also a fluorescence indicator used in experiments acts 
as an exogenous buffer [52]. In the rest state of the neuron the amount of free 
buffers is still large enough such that during stimulation the dominant part 
of incoming free Ca^+ ions bind to these buffers [17]. The steadily ongoing 
binding and dissociation process of Ca^"'" and buffers strongly influence the 
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resulting free Ca^+ dynamics. This applies not only to the phase of growing 
Ca^+ concentration during stimulation but also to the return of the Ca^+ 
concentration to the base level. Therefore, we include both types of buffers 
explicitly in the model using the following kinetic equations 

^^en,ex ~^en,ex ^en,ex + ^en,ex ^ (^en,ex ~ ^en,ex) (6) 



for the endogenous and exogenous buffer [8,13,53]. Here, 6en and 6ex are the 
concentrations of bound buffers, \P^^ and rig^ are the total concentrations of the 
intracellular buffer proteins. A'c'ji cx ^iid ^i^^ the rate constants in units of 
— and — , respectively. As we do not take diffusion into account, we neglect 
all spatial effects of mobile buffers. Also all effects related to different types 
of buffers are neglected. Note, that the source and the sink terms for the 
buffers (Eq. (6)) are the sink and the source terms, respectively, for the Ca^+ 
dynamics (Eq. (1)) 

-Sen,ex (c) ^^en,ex ■ (7) 



The rate constants are in general not available from experiment. We use the 
realistic assumption that the buffer dynamics takes places on short time scales 
compared to the Ca^+ dynamics, i.e. that the buffer dynamics is always in a 
quasi steady state (rapid buffer approximation [15,17,54]). Note, that we ne- 
glect low affinity endogenous buffers in the present model. The rapid buffer 
approximation has two important consequences. First, on the short time scale 
the boimd endogenous and exogenous buffer concentrations 6en,ex become con- 
stant and can be written as 

c 6° 

feen.ex (c) = '"'^ • (8) 

-'^en,ex ~r C 



This implies that the bound buffer concentration is adapted instantaneously 
to the free Ca^+ concentration at each time point t. Therefore, we have 

^eTex (C) = -|&en,ex (c) = -^6en,ex (c) . (9) 



With the help of Eq. (8) and Eq. (9) the Ca^+ dynamics (Eq. (1)) gets a more 
compact form: 

|c^i^W-/. + L}^^^. (10) 



Ten and Tlx are dimensionless variables: 
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^en (c) 



rex(c) 



with Ken = rf 



with Kp 



k+x 



and is a correction factor in Eq. (10) which depends only on the free 

Ca^"*" concentration. The great technical advantage of this approximation [53] 
is that Eq. (10) only involves the dissociation constants /i'en.ex of the buffers 
instead of both rate constants /cjngx ^enex- general the dissociation 
constants are available from experiment. 



2.3 Ca2+ efflux through PMC A and NCX 



The PMCA and the NCX are proteins which actively transport Ca^+ ions 
through the membrane [55,34] and therefore have to be described differently 
from pores. The PMCA pumps Ca^"*" ions against the electrochemical gradient 
using the energy of ATP molecules. The NCX transports Ca^"*" ions out of 
the neuron by exchanging them with Na"^ ions. The Ca^"*" dependence of the 
PMCA-kinetics observed in experiments is well described by the Hill equation 
[56,57,58]. The activity of the PMCA and the NCX is limited to a maximum 
rate [59,34,35] which follows from the molecular structure of the proteins. 
Correspondingly, the free Ca^"*" efflux in Eq. (1) gets the form 

Je (c) = Pp Qp (c) 7p + px S'x (c) 7x . (12) 



This is again formulated on the single protein level, where pp,x represent the 
specific surface densities of the membrane proteins and /p^x the universal max- 
imum activity rates of the PMCA and the NCX, respectively. For the universal 
activity characteristics Qp^^ we assume Hill equations 

^''^''^^W^HHpp ^"""^ ^''^ = (^tTo^ ^^^^ 



which directly depend on the free Ca^"*" concentration c [60]. ifp,x are the half 
activation concentrations and rip x the Hill coefficients. This is justified in the 
case of PMCA [56] while no explicit measurements exist in the case of NCX. 
Nevertheless, it is a reasonable assumption that the principal behavior is also 
well described by a Hill equation. 
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1 

value 




r 

reiereii(;e8 


VDCC 






5v = 14pS 




[44] 


= -4mV 




[44] 


n = 6.3 mV 




[44] 


r = 1 ms 




[4.3] 


V'eq = i7niV 




[44] 


PMCA 






Tp = 0.27 X 10" 


-20 C 
ms 


[56], [61] 


np = 2 




[56] 


Hp = 0.09 ijM 




[56] 


NCX 






7x = 0.48 X 10- 


-18^ 
ms 


[62] 


nx = 1 




see text 


i?x = 1-8 /xM 




[34] 



Table 1 

The values of universal model parameters. 
2.4 Universal properties of the model 

The classification of the model parameters represents an important element 
of our model because it enables us to separate those properties which are 
basically independent of the cell-type from those vi^hich are sensitive for differ- 
ences betM^een cell-types, compartments, and experimental setups. The univer- 
sal properties of the proteins are postulated to remain invariant for all kinds of 
neurons and extracted these values from various single-protein experimen- 
tal results. The values of universal parameters for VDCC, PMCA and NCX 
as used in the model are summarized in Tab. 1. 

VDCC. These channel-proteins are multi-subunit complexes which form a 
voltage sensitive transmembrane pore. Six types of VDCC are known, cor- 
responding to the pharmacological properties (i.e. L-, N-, P/Q-, R-type as 
high-voltage-activated (HVA) channels and T-type as low-voltage-activated 
(LVA) channels. These types are primarily characterized by genes encoding a 
different subunit (e.g. the ai subunit) but an overall matching percentage of 
nearly 50% (for HVA channels only) [63]. The Ca^+ influx into the presynap- 
tic terminal is dominated by the P/Q- and N-type channels [23]. Each type 
of VDCC should have the same biophysical properties in all kinds of tissues 
and animals. Unfortunately the actual precision of the experiments, especially 
taking into account the variety of experimental setups, does not allow to quan- 
titatively identify the biophysical properties corresponding to each channel 
type [64,65,43]. The resulting uncertainties dominate the differences between 
the channel types. Therefore, we chose the values from one single channel 
measurement of N-type channels considering it as a representative member 
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for the VDCC. However, if more conclusive values for different VDCC-types 
are found in future we may disentangle the influence of different sub-types 
on Ca^"*" transients. Concerning the measurements of the P/Q- and N-type 
channels two additional problems exist. Most measurements were performed 
with Ba^+ as carrier ion and the obtained conductivity differs from the Ca^"^ 
conductivity [64]. Additionally all measurements were made at room temper- 
ature (20 — 25°C) and the biophysical properties are temperature dependent 
[66,67]. Thus, the obtained Ca-^+ influx through the VDCC differs between 
measurements made at room temperature from those which were made at 
blood temperature. We neglect both deviations because we could neither find 
convincing extrapolations of the temperature behaviour nor a sufficient rela- 
tion of diff'erent charge carrier ions. Anyhow, those deviations stay within the 
current experimental accuracy. 

PMCA. There exist four different genes which encode the PMCA (PMCAl- 
4) and all of them occur in different splicing variants (labeled by small letters) . 
These iso-forms have different kinetic properties but again we assume that one 
iso-form should have the same properties in each type of tissues and animal. 
The dominant iso-forms of PMCA in the rat brain are the la-, 2a-, 2b-, 3a- 
, 3b-, and 4b-type [32,68,69]. Similarly to the VDCC the measured kinetic 
properties vary considerably due to the different experimental conditions for 
one type of PMCA [70,56,71] and we chose as representative member the 
PMCA2a protein. The Hill coefficient and the concentration of half activation 
are taken from a measurement for rat PMCA [56] and the maximum pump 
current is calculated from [56,61]. We neglect other regulation mechanisms of 
the PMCA-activity e.g. the roles of calmodulin or of ATP (which is considered 
to be available to a sufficient degree). 

NCX. In the case of the NCX three different genes (NCXl-3) with different 
splicing variants are known. Here, we do not specify the type of NCX. We 
chose average values for the maximum current [62] and for the concentration 
of half activation [34], and set the Hill coefficient of the NCX to one. It is 
worth emphasizing that our results are not significantly altered for other Hill 
coefficient. Other regulation mechanisms (e.g. the dependence on Na^) are 
not considered in the model. 



2. 5 Specific properties of the model 

Per definition the specific parameters have to be adjusted to the specific neu- 
ron type and compartment geometry in the experiment under consideration. 
Especially, this concerns the rest state properties (i.e. geometry factor G, rest 
state membrane potential Vq, rest state intracellular Cq and extracellular Cext 
Ca^"^ concentration) which in general are indicated in the experiments. This 
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applies not necessarily to the dissociation constant and the total concen- 
tration mo of the endogenous buffer. Here we chose Calmodulin as represen- 
tative member of buffer molecules and we assume that the four buffer binding 
sites are independent. Our aim is to quantitatively compare the results of our 
model with experiments. Because all visualizations of Ca^"^ dynamics in neural 
compartments are performed with fluorescence indicators we have to include 
the parameters specific for the used type of indicator. The corresponding dis- 
sociation constant can be extracted from corresponding data sheets and the 
indicator concentration is usually indicated in the experimental protocol. 

The most critical specific parameters are the membrane protein surface densi- 
ties (e.g. VDCC Pv: PMCA pp). They are rarely known for specific cell types. 
Even more seldom they are indicated in specific experiments. In addition mea- 
surements of protein densities are usually valid for very specific areas on the 
cell membrane. It is difficult to extrapolate those densities to an average den- 
sity on a whole presynapse as it is used in the model. As our main interest 
points towards the analysis of Ca^+ transients in response to various stimu- 
lation protocols we proceed as follows. We use single action potential Ca^"*" 
transients measured as fiuorescence signal to fit the unknown protein densi- 
ties. Then the values found are retained and used for the analysis of Ca^"*" 
transients after other stimulation protocols. This procedure is explained in 
more details in section 3.1. 



2. 6 Equilibrium, stability and stimulation 

To obtain a realistic simulation of Ca^+ dynamics the model has to recover 
a stable equilibrium state after some stimulation. That is achieved by the 
requirement that in the rest state of the neuron (defined by the rest membrane 
potential Vq and the free rest Ca^"*" concentration cq) the leak conductivity of 
the membrane L exactly compensates the netto current densities coming from 
VDCC (Ji), PMCA and NCX (Je). 

After application of any type of stimulation the free Ca^"*" concentration indeed 
returns to the rest concentration. We did not find any instabilities of the Ca^"*" 
dynamics within physiologically relevant stimulation protocols. The model 
results are robust against variation of any parameters on the qualitative level. 
Quantitatively the results are most sensitive to variations of the membrane 
protein surface densities and of the dissociation constant of the endogenous 
buffer. 

For the stimulation protocols (single action potential and trains of action po- 
tentials) we use membrane potentials V (t) which are simulated with a system 
of coupled differential equations (Hodgkin-Huxley like, not shown here). The 
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parameters of those equations have to be adjusted such that the stimulation 
shapes as found in the experiment under consideration are appropriately re- 
produced. 



3 Results 

To check if our model correctly describes the presynaptic Ca^+ dynamics in 
response to single action potentials, we tuned the model with respect to an 
experiment measuring the Ca^"*" dynamics in presynaptic boutons of pyramidal 
neurons in the neocortex of rats [23]. The corresponding Ca^"*" transients can 
be reproduced for reasonable membrane protein surface densities. Using those 
parameters we can calculate the Ca^+ dynamics due to a 10 Hz tetanus and 
compare the result to corresponding measurements [23] (see section 3.1). A 
qualitative discussion of Ca^+ transients in response to various stimulation 
protocols is provided in section 3.2. Finally we will quantitatively illustrate 
the influence of fluorescence indicators on Ca^+ dynamics and discuss the 
implications for experiments using fluorescence indicators to measure Ca^+ 
transients (see section 3.3). 

3.1 Quantitative model test 

At first we verify that our model can describe Ca^+ transients evoked by sin- 
gle action potentials in single presynaptic boutons from pyramidal neurons in 
the neocortex of rats on the basis of a very detailed corresponding fluores- 
cence measurements [23]. The universal parameters should especially apply 
for this system and are used as shown in Tab. 1. The specific parameters are 
determined directly and indirectly by the measurement [23] and by using inde- 
pendent experimental results. We calculated the geometry factor with the help 
of a fluorescence figure of the whole bouton (Fig. 3A in [23]). We determined 
the dissociation constant of the endogenous buffer using the observation 
that approximately only 1% of the total incoming intracellular Ca^"*" remains 
free [23]. The surface density of the NCX px was determined with the help 
of a measurement which investigates its ratio to the PMC A density Pp [62]. 
The values for the specific parameters are summarized in Tab. 2. They have 
to be carefully interpreted as an approximation for the neuron type under 
consideration. Please note that only two parameters remained free to adjust 
the model to the experimental data: The densities of the PMCA and of the 
VDCC. Note, that any dependence of channel densities on intracellular Ca^"*" 
[72] has not been considered here. 

In general Ca^"*" measurements in neurons use fiuorescence indicators to vi- 
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value 


references 


geometry (steady state) 




G = 3/0.5 Atm 


[23] 


Vb = -70 mV 


[23] 


Co = 0.1 /iM 


[47] 


Cext = 1.5 mM 


[47] 


protein densities (single action potential) 


/9v = 3.1//xm^ 


adjusted 


/9p = 9200/^m2 


adjusted 


Px = 0.033 * /9p 


[62] 


protein densities (10 Hz tetanus) 




= 3.7///m^ 


adjusted 


= 8300/;um2 


adjusted 


= 0.033 * 


[62] 


endogenous buffer 




i^en = 0.5 


[23] 


ftO^ = 4 * 30 /iM 


[13],[73] 


fluorescence indicator 






[23] 


6°^ = 100, 500 /xM 


[23] 


(A///)Z"^ = 1.5 


[23] 


i^f/ff^T = 2.3 


adjusted 



Table 2 

The values of specific model parameter. 

sualize the Ca^"*" concentration and distribution. As already pointed out such 
fluorescence indicators act as an exogenous bufl^er system which, therefore, has 
been included into the model prescription (see section 2.2). With the help of a 
1 : 1 complexation ansatz [23,74,75] we translated the simulation of the Ca^+ 
dynamics into a simulation of the fluorescence signal 

A/// = (A///Ui|;^. (14) 

A/// is the relative fluorescence change, (A///)^^^^ represents its maximum, 
and is the dissociation constant of the indicator. The properties of the 
used indicators (Magnesium Green (Molecular Probes)) are shown in Tab. 2. 
Note, that in the experiment different concentrations of indicator have been 
used [23]. 

The only remaining free parameters are the surface densities of the VDCC Pv 
and the PMC A Pp. We simulate the amplitude and the shape of the action 
potential as applied to the presynapse in the experiment and fit the protein 
densities (see Tab. 2) such that the simulation result agrees with the measured 



13 



80 
60 
40 
> 20 


& -20 



-50 50 100 150 200 250 

0,2 

0,15 

ti; 0,1 

^0,05 





"""-50 50 100 150 200 250 

t [ms] 

Fig. 1. The Ca^+ transient (shown as relative fluorescence change with 100 //M Mag- 
nesium Green) evoked by single action potentials in single boutons of pyramidal 
neurons in the rat neocortex. The upper panel shows the simulated action poten- 
tial compared to the one used in experiment (stars). The lower panel compare the 
experimental results (with kind permission of H.J. Koester and B. Sakmann [23], 
Fig. 15A) to the model results. The parameters as listed in tables Tab. 1 and Tab. 2 
have been used. 
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fluorescence signal evoked by a single action potential (Fig. 15C in [23]). Note 
that the experiment has been executed with 100 iiM Magnesium Green and 
(A///)^°^^'^ has to be chosen correspondingly (see Tab. 2). The result is 
shown in Fig. 1. The amplitude and the shape of the simulated fluorescence 
response are in good agreement with the data observed in experiment. 

The next step is to investigate whether the same model without changing the 
model assumptions is able to describe Ca^^ transients evoked by more complex 
stimulation protocols. We compare the fluorescence response evoked by a 10 Hz 
tetanus to the model predictions. To this end the action potentials in the model 
have to be adjusted to those used in experiment (Fig. 9; upper panel in [23]). 
The experiment has been performed with 500 /iM Magnesium Green and we 
corrected the indicator concentration correspondingly. The maximal relative 
fluorescence change (A///)^°^^^ has not been stated explicitly. Therefore, 
we determine this value by the requirement that the Ca^"*" signal in response 
to a single action potential (Fig. 15C; middle column in [23]) is reproduced 
correctly (data not showai). Finally, the protein surface densities are adjusted 
to reproduce the Ca^"*" response to the first action potential of the tetanus (see 
Tab. 2). The subsequent Ca^"*" signal as predicted by the model turns out to 
be in perfect agreement with the fluorescence signal (see Fig. 2). Note, that 
all others parameters of the model remained unchanged. We conclude that the 
model once fitted to the Ca^+ signal in response to single action potentials in 
a specific type of neuron describes the Ca^+ signal due to tetanus stimulation 
without further changes. 
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Fig. 2. The Ca^"*" transients (shown as relative fluorescence change with 500 /7,M 
Magnesium Green) evoked by a 10 Hz tetanus in single boutons of pyramidal neu- 
rons in the rat neocortex. The measured Ca^+ signal observed in experiment (with 
kind permission of H. J. Koester and B. Sakmann [23], Fig. 9) is compared to the 
corresponding model results (lower panel). The upper panel shows the simulated 
stimulation protocol. The parameters as listed in Tab. 1 and Tab. 2 have been used. 
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Fig. 3. The Ca^"*" transients (without indicator) evoked by a 2 Hz, 3.5s tetanus. All 
other parameters as in Tab. 1 and Tab. 2. The Ca^^ response to each action potential 
stimulation remain independent from each other. 



3.2 Ca^''" dynamics in response to different stimulation protocols 



The main question we address is whether there are characteristic differences 
between the different Ca^"*" transients or not. In the low frequency domain 
(2 Hz, 3.5 s) the Ca^+ response appears as a train of independent single action 
potential responses (Fig. 3). The Ca^+ concentration follows the activation by 
each action potential. Shape and amplitude of each Ca^+ response is unaltered 
compared to the response to a single action potential. This frequencies range 
is typically used for LTD induction [76]. 
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Fig. 4. The Ca^"*" transients (without indicator) evoked by a 20 Hz, 350 ms tetanus. 
All other parameters as in Tab. 1 and Tab. 2. A new baseline of the Ca^"*" concen- 
tration arises during stimulation. 

For higher frequencies (20 Hz, 350 ms) a new baseline in the Ca^"'' concentra- 
tion arises during stimulation (Fig. 4). The Ca^"*" pumps and the Na^/Ca^"*" 
exchangers have not enough time to return the neuron into its rest state. This 
implies an overlap of Ca^+ response to subsequent action potentials. Note, 
that this qualitative behavior has been observed for example in experiments 
performed on dendritic spines of pyramidal neurons [52]. However, we can not 
expect that the presynaptic model also quantitatively describes the Ca^+ re- 
sponse in dendritic spines correctly (especially without including diffusion in 
the model). 

This strong increase of the peak and baseline and especially of the average 
Ca^"*" concentration for stimulations with higher frequencies is a significant 
difference compared to the Ca^+ response to a low-frequency tetanus. One 
may suspect that this behaviour is an important part of a process to change 
synaptic efficiency. For example an accumulation of residual free Ca^"'' is re- 
sponsible for synaptic facilitation in presynaptic boutons [9,77]. Also for the 
induction of LTP it is well established fact, that the intensity of stimulation 
protocols have to overcome some threshold [76]. The strong increase of the 
Ca^+ concentration emerges for frequencies that are in the same range as 
this threshold [78,79]. This hypothesis also agrees with the experimental fact 
that a modest Ca^"^ influx causes LTD, whereas a large Ca^"*" influx triggers 
LTP [76,80]. In general the role of enhanced average Ca^^ concentrations for 
changes of synaptic efficiency has been pointed out by [81]. In order to look for 
a specific key of a concrete induction process of one of the different forms of 
synaptic plasticity one has to consider the synapse in more detail. For instance 
exocytosis steps [82], buffer properties [83], diffusion effects [84] and protein 
kineses [85] are important elements of the induction machinery. This is a task 
for future modeling work. 

During the application of a high-frequency tetanus (50 Hz, 1 s) the induced 
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Fig. 5. The Ca^+ transients (without indicator) evoked by a 50 Hz, Is tetanus. All 
other parameters as in Tab. 1 and Tab. 2. There is a more pronounced enhancement 
of the Ca^"'" baseline. 

Ca^"^ spikes remain in phase with the stimulating membrane potential. The 
enhancement of the Ca^"*" baseline becomes even more pronounced (Fig. 5). 
Again the new baseline saturates during stimulation but on a higher level 
compared to the 20 Hz stimulation. The exact shape of the Ca^"*" signal is a 
result of an interplay of membrane proteins and Ca^+ buffers. However, we 
observed some relations between the general appearance of the Ca^"'' signal 
and some specific neuron properties. The surface density of VDCC basically 
determines the amplitude of the Ca^"*" spike on the top of the baseline. The 
dissociation constant determines the height of the new baseline. The time 
scale on which the Ca^+ concentration returns to its rest state is governed 
by the surface densities of the PMCA, NCX and the dissociation constant of 
the buffer. We emphasize that even higher stimulation frequencies lead to still 
higher baseline levels of the Ca^"*" concentration. The reached baseline grows 
quadratically with the stimulation frequency up to 100 Hz. 

3.3 Fluorescence indicators disturb the intrinsic Ca^"*" dynamics 

In section 3.2 we have seen that for stimulations with frequencies above some 
threshold frequency a new baseline of the Ca^+ concentration arises. In the 
following we show that the use of fluorescence indicators for the visualization 
of Ca^+ in experiments alters the Ca^+ response and especially the emergence 
of a new baseline considerably. In Fig. 6 the Ca^"*" transients evoked by single 
action potentials with and without fluorescence indicator is compared (100 /xM 
Magnesium Green). The amplitude of the Ca^"*" signal is decreased if the flu- 
orescence indicator is used. In addition the Ca^+ relaxation time (needed to 
recover the rest Ca^+ concentration) becomes considerably larger using the flu- 
orescence indicator. These relations are not surprising as the indicator binds 
an important part of incoming free Ca^"^ ions and reduces the free Ca^'*' con- 
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Fig. 6. Ca^"*" transients evoked by a single action potential. The upper panel shows 
the Ca^"*" transient for a neuron with 100 /xM Magnesium Green, and the lower panel 
the Ca^"*" transient which would occur without indicator. The indicator changes the 
amplitude and the shape of the Ca^"*" transient considerably. 



centration. Note, that the larger relaxation times are not caused by a decreased 
pump activity which is induced by the lower Ca^"*" peak. This has consequences 
for the interpretation of experiments using fluorescence indicators. 



The fact that fluorescence indicators act as an additional intracellular buffer 
changing the shape of the calcium transient has been pointed out before 
[15,16]. In Fig. 7 the influence of the fluorescence indicator (500 ^uM Magne- 
sium Green) on the Ca^+ transients evoked by a 10 Hz tetanus in a pyramidal 
neuron of the rat neocortex [23] is shown. The left panel (with fluorescence 
indicator) clearly exhibits a new baseline of the Ca^"^ concentration whereas 
the baseline remains practically unchanged in the right panel (without fluo- 
rescence indicator). Because of the longer relaxation time for the Ca^+ spikes 
the fluorescence indicator facilitates the emergence of a new baseline. Using 
a fluorescence indicator the Ca^"^ spikes overlap already for lower frequencies. 
Consequently, the threshold stimulation frequency for the emergence of a new 
baseline is higher for neurons without indicator compared to neurons that has 
been treated with indicator. A interesting conclusion from that fact is that 
the threshold frequency for the appearance of enhanced Ca^+ level changes 
to lower frequencies if a buffer is present. This may have implications for the 
induction of LTP in experiments using fluorescence indicators. Note, that the 
right panel in Fig. 7 provides a quantitative prediction of the Ca^+ signal 
in the pyramidal neuron of the neocortex as it would have been without the 
use of Ca^+ indicator starting from the Ca^"*" measurement [23]. The same 
reconstruction procedure may be applied to other experiments in future. 
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Fig. 7. Ca^'^ transients evoked by a 10 Hz tetanus. The left panel shows the Ca^+ 
transient with 500 /liM Magnesium Green. The right panel shows the Ca^+ transient 
without indicator. All other parameters as in Tab. 1 and Tab. 2. The presence 
of an indicator changes the Ca^+ signal and the height of the emerging baseline 
considerably. 

4 Discussion 



We developed a model for the presynaptic Ca^"*" dynamics which includes 
VDCC, PMCA, and NCX, as well as endogenous buffers and fluorescence 
indicators. Those ingredients turned out to be sufficient to quantitatively re- 
produce a fluorescence measurement of intracellular free Ca^+ transients in 
response to single action potentials in pyramidal neurons of the rat neocortex. 
On the basis of the parameter set determined by the single action potential 
Ca^"^ transients the model predicts the Ca^"*" response to 10 Hz tetanus (using 
the same neuron type) in a quantitatively accurate way. We conclude that 
from the point of view of the model Ca^+ induced Ca^+-release is not neces- 
sarily involved into the presynaptic Ca^+ dynamics at least for this specific 
type of neuron. 

We would like to emphasize that the single action potential Ca^+ response 
has been produced by fitting only two parameters, i.e. the average membrane 
protein surface densities of the VDCC and the PMCA. All other parameters 
have been determined either by the experiment itself or by independent exper- 
iments. The resulting values for the VDCC density are in good agreement with 
experimental measured values [43]. The PMCA density in contrast turn out to 
be rather high but remain within a range of densities that has been observed 
in experiment [86]. One may think of inverting the line of argumentation and 
interprete the resulting protein densities as prediction for an average density 
on the whole bouton. In view of the difficulties to measure such densities this 
provides an attractive possibility. 

The classification of parameters into universal and specific ones opens the pos- 
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sibility to adjust the model to other experiments without altering the model 
pillars. The (universal) properties of the VDCC for example have been de- 
termined using corresponding single channel experiments. In this context one 
may think of the necessity to include more than one HVA channel type into 
the model instead of restricting oneself to one representative type. However, 
a corresponding analysis revealed that the resulting Ca^+ transients remain 
basically unchanged using two different HVA channels. We conclude that the 
measurement of Ca^"*" transients does not allow to distinguish between differ- 
ent HVA channel subtypes. This may be different if considering LVA channels 
in the model. 

We qualitatively investigated the general behaviour of Ca^+ transients in re- 
sponse to different stimulation protocols. For 2 Hz tetanus we found no char- 
acteristic feature of the Ca^+ signal. This agrees with the interpretation, that 
synaptic facilitation is triggered by calcium enhancements, while synaptic 
depression is a consequence of reduced release site activity [87]. We found 
a characteristic difference between Ca^+ transients in response to 2 Hz and 
50 Hz stimulations. At some threshold stimulation frequency the Ca^+ signal 
does not return to its rest state level after each spike. Instead a new baseline 
emerges on a higher level and the Ca^"^ spikes develop on top of this baseline. 
This general behaviour has also been found in postsynaptic experiments be- 
fore [52]. Indeed, the strong increase in the Ca^+ signal for above threshold 
stimulation frequencies can be interpreted to correspond to the stimulation 
threshold for changing synaptic efficiency [10]. In other words, the emergence 
of an enhanced baseline and, as a consequence, of an enhanced average value 
and increased peak values of the Ca^^ concentration may be thought of as a 
necessary requirement for induction of changes in the synaptic efficiency. The 
height of the baseline increases quadratically with increasing stimulation fre- 
quencies up to frequencies for which the buffer becomes saturated with Ca^^. 
For higher frequencies the baseline increases only linearly. 

Already on the level of calcium transients averaged over the whole bouton 
we can confirm the result, that the endogenous and exogenous buffer concen- 
trations as well as their saturation properties crucially shift the threshold of 
the stimulation frequency [10,15,16]. Especially, it is important to realize the 
strong influence of indicators on the intracellular free Ca^"*" signal [88,89,90]. 
Therefore, we investigated this problem in the framework of the model in 
more details. On the level of single action potentials the use of indicator leads 
to considerably smaller Ca^+ spikes which relaxed more slowly to the rest 
concentration. This implies that the use of indicator shifts the above men- 
tioned threshold stimulation frequency. Therefore, a quantitative evaluation 
of fluorescence measurements should always include the effect of the indicator 
[91,92]. The model presented here, provides a tool to reconstruct the aver- 
aged Ca^+ transient as it would have bean without indicator in fluorescence 
measurements. This has been done in this article for the experiment of [23]. 
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It would be interesting in future to make use of the parameter classification, 
to enlarge the model to calcium and buffer diffusion, and to apply the model 
to other experiments, especially when using other neuron types. A compar- 
ative analysis of presynaptic Ca^+ transients in different neurons may reveal 
neuron-type specific characteristics of the Ca^"'" dynamics. 
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